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Abstract. We present results of numerical simulations of wind variability in 
the magnetic Bl IVe star /3 Cephei. 2D-MHD simulations are used to determine 
the structure of the wind. From these wind models we calculate line profiles for 
different aspect angles to simulate rotation. The results are compared with the 
observed UV wind line profiles. 



1. Introduction 

Winds of many early-type stars show cyclic, or even strictly periodic, variability 
on a rotational timescale. This behaviour was extensively studied in the UV with 
the lUE satellite. Both pulsations and magnetic fields have been suggested as 
the cause for this variability. It is clear that kG magnetic fields as observed in the 
Bp stars can explain the observed, strictly periodic, variability. For other early- 
type stars where no such fields have been detected and which can have much 
denser winds, the origin of the variability is still unknown. Recently, magnetic 
fields have been discovered in several early B type stars, with Cep among 
them . These stars have we aker fields (a few hundred Gauss - iHenrichs et al.l 
I2OOOI : iNeiner et al.l l2003abch and stronger winds than the magnetically strong 
Bp stars, which are mostly of a later spectral type. Using numerical simulations, 
we investigate the effects of the magnetic field on the stellar wind of (3 Cep. 



2. Method 

Our method consists of two steps. First we compute the wind geometry in the 
presence of a magnetic field, using a MHD code that starts out from the relaxed 
non-magnetic wind solution and a dipole field and is evolved in time until quasi- 
steady behaviour is observed. Second, from the end state of our simulations 
we calculate line profiles using a SEI (Sobolev with Exact Integration) method. 
Stellar rotation is simulated b y adjusting the line-of-sigh t of the observer. 

The 2D-MHD code fsee lud-Doula fc Owockill2002l l is based on ZEUS. It 
uses the standard MHD equations, the energy equation with a cooling term ap- 
propriate for stellar winds, and radial CAK-based radiative driving. Since 
Cep is a slow rotator, rotation can be ignored. For our models of Cep, we use 
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Table 1. Stellar parameters of /3 Cep 



spectral type 


BllVe 


polar field 


360(30) G 


mass 


12 Mo 


rotation period 


12.00 days 


radius 


6.5 R© 


inclination 


60°- 90° 


temperature 


-26000 K 




90° 



Figure 1. Results of 2D-MHD simulations of the stellar wind in the presence 
of a 360 G dipole field and a mass- loss rate of the order of 10~^ M© yr~^. 
One can see the outfiowing "disk" , in the magnetic equator, which has higher 
temperatures, higher densities and lower velocities than its surroundings. 
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basic stellar parameters as listed in Tabled We treat mass loss rate, M as a 
parameter that ranges from ^10~^-10~^ Mq yr~^. The line profiles are calcu- 
lated using the SEI method. Opacities are determined using solar abundances 
and a (recently included) LTE ionisation balance. At low temperatures (^ Teff ) 
the CIV density is increased with a few percent due to the Auger effect. The 
highest relative CIV densities are found in shock heated sheets of ~10^ K. Our 
calculations do not account for the combined effect of the doublet members. 



3. Results and Conclusions 

For a model with a polar field (Bp) of 360 G and mass loss rate of ^10~^ 
Mq yr~^, Fig. d shows the density, temperature, and velocity at a typical time 
snapshot long after any influence of the initial state has faded. The material is 
guided by the magnetic field towards the magnetic equator, where it is heated 
by shocks. Due to the higher densities in this "disk", the radiation force is 
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Figure 2. (Left) Strictly periodic UV wind-line variability as observed in 
the CIV doublet of P Cep. The top panel shows the significance of the vari- 
ability. (Middle) Results for a model with Bp =190 G and a fixed CIV abun- 
dance. The basic features are reproduced, but fine-tuning of the parameters is 
needed. (Right) Results for Bp=360 G, including a temperature dependent 
CIV abundance. Note the 0.25 offset in phase relative to the observations. 



HD 20502 l/p Cep CIV 







reduced and radial velocities are relatively small. In Fig.I21we show the observed 
CIV variability as a function of the rotational phase and calculations for the 
blue line of the doublet. Both in the observations and the simulations with a 
fixed abundance, most absorption is seen along the magnetic equatorial plane. 
For the simulations including the ionisation balance, the pattern appears to be 
offset by 0.25 in phase. This is caused by the high temperature in the disk which 
reduces the CIV abundance. A more detailed treatment of the Auger effect is 
likely required to account for this offset, since most X-rays will come from the 
shock heated disk. Recent results (Leone et al] l2005r ) show that magnetic fields 
measured with the LSD method could be overestimated. If f3 Cep would indeed 
have a weaker field, this could help to better explain the observed wind terminal 
speed. Although more work needs to be done, in principle the effect of the 
magnetic field on the wind of /? Cep can explain the observed UV line-profile 
variability. Further modelling of other magnetic early B stars and wind lines 
like NV and SilV, will allow a detailed view on how stellar winds are affected 
by magnetic fields, and the importance of the Auger effect. 
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